Critical temperature and giant isotope effect in presence of paramagnons. 
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We reconsider the long-standing problem of the effect of spin fluctuations on the critical tem- 
perature and isotope effect in a phonon-mediated superconductor. Although the general physics 
of the interplay between phonons and paramagnons had been rather well understood, the existing 
approximate formulas fail to describe the correct behavior of T c for general phonon and paramagnon 
spectra. Using a controllable approximation, we derive an analytical formula for T c which agrees 
well with exact numerical solutions of the Eliashberg equations for a broad range of parameters. 
Based on both numerical and analytical results, we predict a strong enhancement of the isotope 
effect when the frequencies of spin fluctuation and phonons are of the same order. This effect may 
have important consequences for near-magnetic superconductors such as MgCNi3. 

PACS numbers: 74.20.Mn, 74.62.-c, 74.70.-b 



In the last decade a large number of superconduc- 
tors were discovered in which enhanced spin fluctuations 
(SF) play a role in the superconductivity, e.g., Sr2Ru04, 
MgCNi3, e-Fe, ZrZn2, and others, bringing about new 
and interesting physics. However, understanding such 
materials, even at an intuitive level, has been hindered 
by the lack of a simple formula that would approximate 
the full Eliashberg theory in a compact analytical form, 
as the conventional McMillan formula (MMF) does. As 
a result, uncritical generalizations of the latter have been 
used as a substitute, despite the fact that, as we will show 
below, some of them are too approximate or outright in- 
correct. In this Letter we present an analogue of the 
MMF, derived in a controllable way and tested against 
numerical solutions of full Eliashberg equations, includ- 
ing interaction with SF (paramagnons). We point out 
the possibility of a giant phonon isotope effect induced 
by SF. We will also apply this theory, as an example, to 
a nearly- ferromagnetic superconductor, MgCNi3. 

The understanding that SF are pair breakers in conven- 
tional superconductors is nearly as old as the BCS theory 
itself 0. Moreover, it was soon realized that strong cou- 
pling manifests itself in a nontrivial way in the presence 
of SF 0, ■ In a number of papers numerical solutions 
of the Eliashberg equations were presented, incorporating 
phonon a 2 F p (uj) as well as SF a 2 F s (u>) spectral functions 
(see, e.g., Refs. |j|). However, solving the full Eliashberg 
equation is not always an option, and does not provide 
as much physical insight as analytical treatment. An an- 
alytical tool comparable to the famed MMF is needed. 

Retrospectively, one can realize that the overwhelm- 
ing success of the MMF is due to three facts: (a) it can 
be derived analytically using simple approximations, (b) 
it includes Coulomb repulsion effects, (c) it has three 
universal adjustable parameters, which, after little tun- 



ing, produce an expression which is surprisingly accu- 
rate for a large range of phonon frequencies and coupling 
strengths. Compared to the BCS equation, the MMF in- 
cludes three essential pieces of additional physics: effec- 
tive mass renormalization, logarithmic reduction of the 
Coulomb repulsion, and proper (logarithmic) averaging 
of the phonon frequency. All three effects can be de- 
rived analytically in some approximations. In fact, it 
is known that the functional form of the MMF can be 
derived in two different ways. One, known as the square- 
well model j5j, uses the Matsubara representation, where 
the coupling with the phonons is parametrized in terms 
of the matrix A(n, n'). The model assumes two different 
approximations for the same function A(n,n'), depend- 
ing on whether it is used in the equation for the mass 
renormalization Z or in the one for the gap function <f>: 

Xz(n,n') = \ p Q(u p - |w„_„/|) (1) 
A^(n, n) = \p@(wp—\u n \)®(wp—\w n >\). 

This models leads to an equation for the critical temper- 
ature, T c , 

T c = aw log exp{-6(l + A Z )/[A - + cX z )}, (2) 

where the theoretical parameters are a = 1.14, b — 
c = 1, Xz = A<£ = X p = 2 f uj~ 1 a 2 F p (uj)duj and 
Aplnwi g = 2 J °° iv^ 1 In oj a 2 F p (uj)duj. The renormal- 
ized Coulomb potential is reduced from its bare value 
fi as fi* = /i/(l + /iln^ 2 -), where u>c characterizes 
the frequency cutoff of the Coulomb interaction. The 
MMF formula is given by Eq|21with optimized parame- 
ters a = 1/1.2, b = 1.04, and c = 0.62. 

SF, as opposed to phonons, induce repulsion for sin- 
glet pairs. However, they contribute to the mass renor- 
malization just the same. Therefore the first instinct is 
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to let Xz = \ p + A s , where A s describes the SF, and 
A</> = A p — A s . Eq|2]with this modification and standard 
a, b and c is the one routinely used in the literature for 
materials with SF (e.g., Refs. [113.11 )• 

Obviously, using two different approximations for the 
same physical function X(n,n') depending on whether it 
appears in the first or second Eliashberg equation can- 
not be justified by any logic. It appears that the MMF 
formula can be fortuitously derived in this way, but, as 
we will see below, this approach fails when SF are in- 
cluded. An alternative derivation of the MMF utilizes 
the real frequency axis formalism 0. The one-mode ap- 
proximation is used, which assumes an Einstein phonon 
at a frequency u> p , i.e., a 2 F(u>) — X p ujpS(uj — lo p )/2. The 
Eliashberg equations are then solved iteratively. After 
the first iteration one obtains 
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T c = I.Mwpexp < — — — 
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which is similar to the square well formula Eq|2with a = 
1.14/ y/e = 1/1.44 (note that this value of a is much closer 
to the optimized one), b = 1, and c = 0.5/ (1 + A p ). This 
approach is a controllable approximation with a concrete 
physical meaning. However, it has never been applied to 
superconductors with SF. 

On the contrary, several attempts to apply the square 
well model to SF have been reported. In Refs. 0,0 the 
following expression was derived (for fi* = 0) : 

T c = l.l^X^ exp{-(l + A p + A s )/(A p -A s )], (4) 



with v = X p /(Xp — A s ) 



A^ 
Ap — A s 



Ad — A s 
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1 + Ap + A a 
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where the choice (a) is due to Carbotte et al. Q, and 
(b) to Vonsovskii et al. 3]. Unfortunately, neither au- 
thors give details of their derivations, so we do not know 
what was different in their models. We were not able to 
reproduce either result. The latest paper utililizing the 
square well model (in the weak coupling limit) is that by 
Shimahara [1(3 ] . Our own result for the square well model 
reduces to that of Ref.IlOl in the weak limit, and reads 



T c = 1.14w p exp[— - 



1 + A„ + A 



p 



T c = l.lAu) a exp[— 



A P 1+A S +A 3 ln^ 
1 + X s + X v 



X p (l+X p ) 
1 + Ap — A p In - 



A., 



LU S > LU p (5) 



U s < LU p (6) 



Unlike EqU Eqsgl El reduce to the McMillan form upon 
substitution u> s — > loc 3> lo p , X s — > fi, as it should. 

Given the controversy about the square-well model, 
it is desirable to have a derivation in a controllable ap- 
proximation, such as the real frequency axis formalism 



of Refs. 0] . Assuming an Einstein phonon at a frequency 
ujp and an "Einstein" paramagnon at io s , 2a 2 F{uo) = 
XpL0 p 5{uj — ujp) — X s lo s 5(lo — co s ), we obtain the following 
iterative solution of the Eliashberg equations. 



T c = l.Ucup' 
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x exp ■ 
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(7) 



In — 



For ujp — > uj s , K = -1/2, and at /i* = 0, Eqdreduces to 
Eq^with v = X P /(X P — X s ). 

As usual, the ultimate test for any approximation is 
numerical calculations. We solved the Eliashberg equa- 
tions for a variety of model a 2 F(ui) including SF and 
compare them with the proposed analytical formulas. In 
Fig Q we show this comparison for the simplest "one- 
mode" approximation, one phonon and one paramagnon 
(we have verified that other model spectra lead to simi- 
lar results). As we can see, while the Eq[7| as well as its 
simplified version Eq0Ja), describe the numerical results 
rather well when u) s and ui p are comparable, the latter 
fails at lu s ujp, and both fail at u s uj p . Both effects 
can be easily understood: Eq^ includes u s in a negative 
power in all regimes, thus leading to a total suppression 
of superconductivity at u> s — > oo. In reality, in this limit 
the negative effect of the SF is renormalized down loga- 
rithmically in the same spirit as the Coulomb repulsion. 

Eqs. 00] diverge at w s -> 0. This is due to the fact 
that the derivations above assume that u) s ,uj p > 7rT c . It 
is possible to treat this regime separately. If uj s <^ T c , 
the SF act as static magnetic defects, and the standard 
theory of the magnetic pair-breaking jllj can be applied. 
In the Matsubara representation, at tv s = one needs 
only to keep the term with n = m in A s (tu n — uJ. Then 
the equations reduce to the standard form 0, with 
the pair-breaking parameter 7 = (1/2tp)/ttT c = X s . In 
the weak coupling limit, T c is 



r c = r e0 exp[V(l/2)-V(l/2 + 7)], 



(8) 



where T c q = T C (X S = 0). One important difference exists 
between pair-breaking by SF with uj s — and by mag- 
netic impurities: in the former case the pair-breaking 
parameter 7 now does not depend on T c . This has con- 
sequences for the isotope effect, as we will see below. 

For small, but finite oj s <ti irT c summation of X s (n — m) 
over n — m provides the expression for the pair breaking 
rate in Eq|S] 7 = A. q ^- coth ^j- . This result coincides 
with Eq. 5.8 of Ref. 12;] for dynamical pair breaking in 
anisotropic superconductors if the anisotropy parameter 
g (as defined in Ref. ^|) i s se t t° "1- When uj s increases, 
T c drops sharply with a complete loss of superconductiv- 
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FIG. 1: Comparison of T c and isotope coefficient with the 
exact numerical calculations, (color on line) 



ity at u s = uj* = e~ c T cQ /2-y (where e c ~ 1.78). How- 
ever, the condition uj s <C 7rT c used in the derivation of 
EqlHl is lost well before w* (in fact, at uj s ~ w*/2). 

One can take into account the strong coupling effects 
in the square-well model, resulting in a renormalization 



7 -> 7/(1 + A P ) 



coth ; 



As the comparison 



1 + A P 2w s 2Ws 

with numerical calculations shows (Fig. this approx- 
imation underestimates T c . However, it illustrates why 
T c flattens out at a finite value smaller than T c n, when 
uj s —> 0, instead of raising as Eq0 suggests. 

We also show in Fig^Jthat both Eq^Jb) and the square 
well model, Eqs[5land|Hl disagree qualitatively with the 
numerical results in the whole range of uj s . 

We will now turn to the isotope effect. Looking at 
EqQJ one observes that the isotope coefficient, f3 = v/2 = 
Xp/2(X P — A s ) > 0.5, is always enhanced compared to 
its BCS value and is independent of the SF frequencies. 
Clearly, this should hold approximately in the range of 



the applicability of this formula, uj s z 
the more accurate Eq[7| yields for (3 



> 7rT c . Indeed, 
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(9) 



F(r) = (r 2 -2rlnr-l)/(r-l) 2 



The second term here is the correction to Eq0] It can be 
of either sign, since with growing r the F(r) monotoni- 
cally grows from -1 to 1, and F(l) — 0. As discussed, Eq0] 
itself becomes invalid at oj s < ttT c . As u> s — > 0, according 
to Eq|Sl j3 — 0.5 (note that in the case of magnetic im- 



purities (3 > 0.5 due to the dependence of 7 on T c 
Therefore, the isotope effect has to have a maximum at 
some < ui s < uip, and /3 max > X p /2(X p — A s ). 

This is confirmed by numerical calculations, which 
do show that the maximum isotope effect for given 
A s , An is achieved close to u> s ~ lu p and is not far from 
Xp/2(X p — A s ). This is a very import result, and we em- 
phasize it again: if superconductivity is depressed by spin- 
fluctuations, the total isotope effect increases compared to 
its BCS value. 

We shall now apply this formalism to a superconduc- 
tor where T c is believed to be substantially suppressed 
by SF, MgCNi 3 HQES, which has attracted substan- 
tial interest not because of its relatively modest criti- 
cal temperature, T c sa 8 K, but because of its unusual 
antiperovskite crystal structure and proximity to ferro- 
magnetic instability. The latter was first pointed out by 
Rosner et al. 0, who believed in such strong coupling 
with SF that they proposed a p— wave superconductivity. 
Singh and Mazin [14J also came to the conclusion that SF 
should play a role in superconductivity of MgCNi 3 , but, 
based on their frozen phonon calculation, they deduced 
a large electron-phonon coupling constant (A p > 1) due 
to the bond-bending Ni phonons. They reconciled this 
relatively large A p with a modest T c within a scenario of 
s-wave phonon-induced superconductivity depressed by 
SF. Later this scenario was re-invented by Shan et al. 
|l5| . who proved the s-symmetry of the order parame- 
ter by tunneling experiments. This point has been since 
confirmed by several groups and seems to be well estab- 
lished. 

Singh and Mazin's 0] prediction of the Ni phonon 
playing the major role in the electron-phonon coupling 
in MgNiC3 was based on a limited number of calcula- 
tions at a high-symmetry point in the Brillouin zone, and 
therefore was more an educated guess than a quantita- 
tive argument. A quantitative analysis was provided by 
Ignatov et al. 0] , who performed linear-response calcu- 
lations of the phonon frequencies and their coupling with 
electrons for the whole Brillouin zone. They found a gi- 
gantic coupling for the Ni bond-bending modes, and the 
most strongly coupled modes (the mode considered by 
Singh and Mazin was not among them) actually unsta- 
ble. In other words, they found a set of double-well type 
instabilities involving mostl y N i atoms. This was verified 
by EXAFS measurements |l6j. Ignatov et al. 0] es- 
timated the total electron-phonon coupling constant as 
1.5 and the logarithmically averaged phonon frequency 
as 131 K. 

Thus, the scenario of Ref. was modified in Ref. 
(lif in the sense that electron-phonon coupling and su- 
perconductivity were coming from highly anharmonic 
predominantly Ni modes, but not exactly the simple ro- 
tations of the Nig octahedra considered in Ref. [lj| . Un- 
fortunately, strong anharmonicity of these modes makes 
it impossible to evaluate their coupling with electrons 
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FIG. 2: T c for the electron-phonon spectral function calcu- 
lated in plU for MgCNi3 (inset), (color on line) 



in the linear response calculations, but it is obviously 
strong. However, one can estimate the electron-phonon 
and electron-paramagnon coupling indirectly from exper- 
imental data. Indeed, specific heat renormalization, from 
different reports, ranges from 2.6 to 3.1 (see Ref. [l?J and 
refs. therein), implying that the sum A p + X s varies be- 
tween 1.6 and 2.1. Walte et al. [ItJ estimated to p s» 143 
K, smaller than, but comparable to the calculation in 
Ref.|l6|. uj s w 25 K and the mass renormalization due 
to paramagnons as 1 + X s » 1.43. Then, using MMF, 
fi* = 0.13, and T c = 6.8 K, as measured for their sam- 
ples, they deduced A p =1.91. 

However, there are several problems with this deriva- 
tion. First of all, as shown above, the proper formula is 
Eq0] Using this formula instead of Eq|3 and keeping all 
their other parameters, we get a much more reasonable 
number, A„=1.61, not far from the value of 1.51 obtained 
in Ref. 0. However, the SF model adopted in Ref. 
cannot be considered as proven. It is based on the dis- 
putable assumption that the upturn of the specific heat 
quotient at low temperature and high magnetic field is 
due to the paramagnon contribution to specific heat, but 
there many other explanations of this effect. 25 K seems 
to be unrealistically soft. Also, low T c and high residual 
resistance cast doubt on the sample quality in this study. 

Here we adopt a different approach: we adopt the cal- 
culated values A p =1.5 and oj p — 131 K, in the har- 
monic approximation, and total mass renormalization 
1 + A p + A s = 2.85, so that A s = 0.35. The results of 
the numerical solution of the Eliashberg equations with 
the a 2 F(ui) function calculated by Ignatov et aL0| and 
fi* = 0.12 are shown in Fig .2, together with the curve 
calculated from Eq[7| This way, we find ui s ~ 50 K, 
which, we believe, is a more realistic number than 25 K. 
The corresponding total isotope effect coefficient is 0.75. 

This may sound in agreement with the recent experi- 
ment by Klimczuk and Cava |l3| , who have measured the 



isotope effect to be 0.54 on carbon only. If the total iso- 
tope effect is 0.75, this suggests a seemingly reasonable Ni 
isotope effect of 0.21, suggesting that Ni phonons couple 
with the electrons twice weaker than C ones. Unfortu- 
nately, the first-principles calculations suggest that the 
Ni modes couple with electrons at least an order of mag- 
nitude stronger than the C modes (there is hardly any 
C character present at the Fermi level). In the moment, 
the only way to reconcile this with the measurements of 
Ref. is to assume that the observed isotope effect is 
not a result of the frequency shift of the C modes, but 
of some subtle changes in the crystal structure induced 
by the isotope substitution. Such a possibility is sug- 
gested by an earlier study 0], where it was found that 
(i) T c depends on the lattice parameter at a rate of « 310 
K/A, which translates an error of ±0.0015 A in the lat- 
tice parameter|l3j into an error of ±0.46 K in T c , larger 
than the isotope shift of 0.3 K, and (ii) that two samples 
with the same lattice parameter and the same neutron- 
measured C content have T c differing by 0.71 K. A possi- 
ble explanation is that, given the proximity of MgCNi3 to 
a ferromagnetic instability, crystallographic defects may 
induce local magnetic moments which, in turn, work as 
pair-breakers. The concentration of such defects, even 
for the same net C content, may depend on the sample 
preparation and, possibly, on isotope substitution. 

Therefore further studies of the isotope effect both on 
C and on Ni are necessary, in particular combined with 
accurate measurements of the isotope shift of the phonon 
modes. 
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